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■ Various approaches to estimate turbulent transport coefficients from numerical sim- 

ulations of hydromagnetic turbulence are discussed. A quantitative comparison between 
the averaged magnetic field obtained from a specific three-dimensional simulation of a ro- 
QQ ' fating turbulent shear flow in a slab and a simple one-dimensional alpha-omega dynamo 

. model is given. A direct determination of transport coefficients is attempted by calculat- 

ing the correlation matrix of different components of the field and its derivatives. This 
matrix relates the electromotive force to physically relevant parameters like the tensor 
components of the alpha-effect and the turbulent diffusivity. The alpha-effect operating 
on the toroidal field is found to be negative and of similar magnitude as the value obtained 
in previous work by correlating the electromotive force with the mean magnetic field. The 
turbulent diffusion of the toroidal field is comparable to the kinematic viscosity that was 
determined earlier by comparing the stress with the shear. However, the turbulent dif- 
^ I fusion of the radial field component is smaller and can even be formally negative. The 

' method is then modified to obtain the spectral dependence of the turbulent transport 

coefficients on the wavenumber. There is evidence for nonlocal behavior in that most 
^ , of the response comes from the smallest wavenumbers corresponding to the largest scale 

' possible in the simulation. Again, the turbulent diffusion coefficient for the radial field 

component is small, or even negative, which is considered unphysical. However, when the 
diffusion tensor is assumed to be diagonal the radial component of the diffusion tensor is 
positive, supporting thus the relevance of a nonlocal approach. Finally, model calculations 
are presented using nonlocal prescriptions of the a-effect and the turbulent diffusion. We 
emphasize that in all cases the electromotive force exhibits a strong stochastic component 
which make the alpha-effect and the turbulent diffusion intrinsically noisy. 



1 Introduction 

In recent years, significant progress has been made in understanding the generation of large 
scale magnetic fields by dynamo action. One of the outstanding questions is whether such large 
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scale dynamos work in the limit of high conductivity and in the presence of finite amplitude 
magnetic fields. These questions have so far mostly been addressed in the framework of fully 
periodic boxes with external forcing (Cattaneo & Vainshtein 1991, Vainshtein & Cattaneo 1992, 
Cattaneo & Hughes 1996, Brandenburg 2001; hereafter referred to as BOl). Subsequently, con- 
cerns have come up, because periodic boxes have the pecuhar property that magnetic helicity is 
conserved. In that case the a-effect may well be 'catastrophically' quenched by finite amplitude 
magnetic fields if the magnetic Reynolds number is large (Blackman & Field 2000, BOl). 

Unfortunately, externally forced simulations with non-periodic, open boundaries do not 
seem to alleviate the catastrophic quenching problem (Brandenburg & Dobler 2001; hereafter 
referred to as BD). This is the reason why there is now an urgent need to consider simulations 
of naturally driven turbulence. Apart from convection, for which the a-effect has been studied 
extensively (Brandenburg et al. 1990, Ossendrijver et al. 2001), shear flow turbulence that is 
driven by the magnetorotational instability is another important example (Brandenburg et al. 
1995, hereafter referred to as BNST95, Hawley et al. 1996). In the case of naturally driven 
turbulence, a-effect and turbulent diffusivity are no longer scalars but tensors, nor are they 
just coefficients but integral kernels. The purpose of this paper is to explore these aspects of 
inhomogeneous turbulence using three-dimensional simulations. 

Simulations relevant to astrophysical bodies or to the Earth that include the effects of rota- 
tion and large scale velocity shear have recently displayed strikingly coherent spatio-temporal 
order (Glatzmaier & Roberts 1995, BNST95). As an example we reproduce in Figure |I| a 
space-time diagram (or butterfly diagram in solar physics) of the mean magnetic field of an 
accretion disc simulation of Brandenburg et al. 1996 (hereafter referred to as BNST96), which 
is an extended run of BNST95. In this particular simulation (Run of BNST96) the symmetry 
of the magnetic field has been restricted to even parity, so the computation has been carried 
out in the upper disc plane, < 2; < L^, where Lz = 2H is the vertical extent of the box, H is 
the gaussian scale height of the hydrostatic equilibrium density, and z = corresponds to the 
equatorial plane. 

It is interesting to note that the spatio-temporal behavior obtained from the three-dimensional 
simulations resembles in many ways what has been obtained earlier using mean-field models. 
In particular, the magnetic fields of the simulations of Glatzmaier & Roberts (1995) show mag- 
netic field patterns with basic dipole symmetry and occasional field reversal, similar to what 
is seen in mean-field dynamo models (e.g. Hollerbach et al. 1992, Jones & Wallace 1992). Of 
course, the simulations of Glatzmaier & Roberts (1995) are much more realistic and show addi- 
tional features such as the magnetic dipole inclination. However, we are not aware of a detailed 
attempt to relate the three-dimensional simulations to mean-field theory. At present there 
are no estimates of the relevant alpha-effect, or even its sign, and whether it is actually the 
alpha-effect that is responsible for the large scale field generation in these simulations. Other 
proposals for causing large-scale field generation include, for example, the incoherent a-effect 
(Vishniac & Brandenburg 1997), negative magnetic diffusivity effects (Zheligovsky, Podvigina, 
& Frisch 2001), or effects that conserve magnetic helicity explicitly (Vishniac & Cho 2001, 
although there is non-supportive evidence; see Arlt & Brandenburg 2001). 

The aim of this paper is to assess the possible relation between simulations and mean-field 
theory using a local simulation of BNST96 that is readily available to us. This simulation has 
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Figure 1: Horizontally averaged radial and toroidal magnetic fields, and By respectively, 
as a function of time and height, as obtained from the fully three-dimensional simulation of 
BNST96. Time is given in units of rotational periods. Trot- (No smoothing in 2; or i is applied.) 
Dotted contours denote negative values. 

the advantage of being in cartesian geometry, and one can define averages that arc naturally 
dependent on just one spatial coordinate (z). Furthermore, geometrical or boundary effects are 
not very pronounced in that model. 

In the particular process considered by BNST95 an initial magnetic field of suitable mag- 
nitude, but vanishing average {B^ = By = Bz = 0), leads after at least thirty rotational 
periods to a self-sustained turbulent state, which is statistically steady and independent of 
the detailed initial conditions. This process results from two instabilities, both with positive 
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feedback. There is the magnetic shearing or Balbus-Hawley (1991) instabihty that generates 
the turbulence and a dynamo instabihty that regenerates the magnetic field. The energy that 
drives the dynamo is constantly being tapped from the large scale velocity shear. 

At first glance it may seem inappropriate to invoke a mean-field description in the present 
case which is highly nonlinear and where the turbulence itself is driven by the dynamo-generated 
magnetic field. However, we adopt here the view that the mechanism that leads to the gener- 
ation of large scale magnetic fields is distinct from the feedback cycle that leads to the smaller 
scale magnetic fields driving the turbulence, which in turn drives the magnetic field. There have 
been some attempts to model this feedback cycle using a set of phenomenological equations 
(Tout & Pringle 1992, Regos 1997), but these models too involve an a-effect that is prescribed 
in an ad hoc fashion. Unfortunately, these models do not reproduce a number of features seen 
in the simulations (ratio of poloidal to toroidal fields, magnitude of the accretion torque, role of 
the Parker instability, and significance of the vertical field.) Also, it should be emphasized that 
the large scale field seen in BNST95 is really a consequence of vertical stratification. Models 
without stratification (Hawley et al. 1996) do not show large scale fields, although the feedback 
cycle of dynamo-generated turbulence still works. This confirms that the large scale field gener- 
ation can indeed be considered as a process operating on top of the otherwise highly nonlinear 
feedback cycle leading to dynamo-generated turbulence. If this is the case one might be able to 
describe the large scale dynamics in terms of averaged equations. This would imply a closure 
relation between the nonlinear induction term (the electromotive force) and the mean magnetic 
field. It is this relation which is at the center of the present paper. 

On the one hand, our findings support basic aspects of the mean-field concept. On the other 
hand, when trying to estimate the mean-field transport coefficients (alpha-effect and turbulent 
diffusion) it becomes necessary to go beyond simple parameterizations, to adopt tensorial and 
nonlocal forms of closures, and to allow for stochastic effects. While these results are primarily 
applicable to local simulations of accretion discs, it will be interesting to make qualitative 
comparisons with simulations of other systems (geodynamo, solar dynamo). 

Our parameterization of transport coefficients depends implicitly on the magnetic field 
strength, because the data used are from a fully nonlinear calculation. However, we are not 
able to investigate this dependence in our present approach, because the time-averaged level of 
magnetic field is rather stable in the model. We should however emphasize that in the present 
model, where the turbulence is magnetically driven, the magnetic energy is strong (three times 
in superequipartition with the turbulent energy) and that the turbulent transport coefficients 
are therefore already strongly affected, and possibly even enhanced, by this strong magnetic 
field. 

Before we begin we summarize a few important properties of the simulation of BNST96. 
The calculation was carried out in a cartesian box where overall rotation, radial linear shear 
and vertical density stratification are included. The boundaries in the toroidal (y) direction are 
periodic and in the radial (x) direction sliding periodic (Hawley et al. 1995), which is periodic 
with respect to positions that shift in time. At the top the boundaries are impenetrable and 
stress free, and the horizontal components of the magnetic field vanish. At the lower boundary 
of the box a symmetry condition is applied. 

We consider the horizontally averaged magnetic field components in the two horizontal 
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directions, and By] see Figure |T]. The averaged vertical field is conserved, because periodic 
(or shearing-periodic) boundary conditions are used in the horizontal directions. Since the 
averaged vertical field vanishes initially, it vanishes at all times, see BNST95. We observe the 
following remarkable properties of the horizontally averaged field: (i) its sign changes every 
approximately 15 rotational periods, (ii) at larger heights the reversals occur somewhat later 
than at lower heights, and (iii) the toroidal and radial fields are out of phase by about Sit /A. 
The latter is more clearly seen in curves showing the averaged field as a function of time (e.g. 
Fig. Ic in BNST96). Property (ii) can be interpreted as a pattern migration away from the 
equatorial plane at speed c > (c ~ 0.024:QH). In addition, there are properties related to 
the symmetry of the field about the equatorial plane; see Brandenburg (1998) for a discussion 
of the dependence on boundary conditions (and the striking agreement with predictions from 
simple afi-dynamos) . 

We begin by discussing the signs and magnitudes of various helicities that are relevant in 
connection with the a-effect. We then compare quantitatively a simple mean field model with 
the result of the simulations. We also consider the direct determination of turbulent transport 
coefficients allowing for the possibility that the transport coefficients may depend either on 
z or on the vertical scale (i.e. the vertical wavenumber). Finally we present a model with a 
scale dependent (or nonlocal) formulation of magnetic diffusion and a-effect, as well as a model 
where strong stochastic fiuctuations are included. 

2 Magnetic and kinetic helicities 

In the light of recent results for a and rjt in forced turbulence, it is useful to discuss some relevant 
properties of the present accretion disc simulations. Firstly, as expected for a rotating stratified 
medium, the kinetic helicity is negative in the upper disc plane. In BOl the magnitude of the 
ratio /Ceff = u) ■ u/u^ was found to be equal to the forcing wavenumber, k{. (Here, a; = V x u 
is vorticity and overbars denote horizontal averaging.) Applying this to the disc simulations 
we find the effective forcing wavenumber between 2 and Ain H < z < 2H, and nearly zero in 
< z < H; see Figure |l|. These wavenumbers correspond to a length scale of about 3, which 
is comparable to the scale of the box. Thus, one cannot expect to see some of the pronounced 
features of BOl associated with an at least modest amount of scale separation. 

In BOl the small scale current helicity, j ■ b, was of the same sign as uj ■ u and of the same 
magnitude. (Here, j = V x 6 is the small scale current density in units where the magnetic 
permeability is unity.) This is not the case in the present simulations. First of all, j ■ b is 
positive near the top of the box in 1.5H < z < 2H, and at most only about 10% of \uj ■ u\. In 
< z < 1.5H the small scale current helicity is essentially fiuctuating about zero. However, 
as in BOl, there is also here a tendency for the total current helicity, J ■ B + j ■ b, to cancel 
to zero, although not quite, and only near the top boundary where large scale and small scale 
current helicities tend to have opposite signs. This is however different from the model studied 
by BD, where open boundary conditions (the same that are used here) were found to yield a 
modest helicity flux out of the domain, which offsets the otherwise strict steady state balance 
between large scale and small scale current helicities that was found in BOl and also here. In 
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Figure 2: Vertical dependence of magnetic helicity H = Hm + 'Hf and magnetic helicity flux 
^ = + ^f, where the subscripts m and f denote the contributions from the mean and 
fluctuating fields. Magnetic helicity is normalized by Tio = J B'^dV, where fci = vr/L^ is 
the smallest wavenumber in the vertical direction. The magnetic helicity flux is normalized by 
J-'o = T-CqQH, where H = Lz/2 is the density scale height of the disc. 



the present simulations the magnetic helicity flux, is however very small; see Figure ^ 

It is somewhat surprising that in the present simulations the resulting magnetic helicity flux 
is so weak. This could be related to the fact that here, in contrast to BD, the turbulence is 
driven self-consistently by the resulting magnetic field. Near the boundaries there are relatively 
sharp gradients in all quantities, which is typical of boundary layer behavior, and is unlikely to 
occur in more realistic cases where free turbulent exchange and flows through the boundaries 
can occur. 

The fact that the small scale current helicity is positive indicates that a is negative. The 
connection between the two is based on a formula given by Keinigs (1983), a = —rjj ■ b/B , 
where rj is the microscopic magnetic diffusivity.^ Another more direct estimate for a comes 
from considering the balance of different induction effects, aB + uxBon the one side and 
iv + Vt)J on the other. Taking the dot product with B shows that a = +rjtJ ■ B/B . Thus, 
since J ■ B is found to be negative we again find that a is negative in the upper disc plane. 
(The negative sign of J ■ S is in agreement with a negative sign of the magnetic helicity Ti 
throughout most of the computational domain; see Figure |^.) 

Perhaps the most convincing explanation for the negative a is that intense parts of a flux 
tube contract (to maintain pressure balance along field lines), but are also most buoyant. If this 
contraction is stronger than the expansion associated with the rise into a less dense medium, 

2 

^We note that the resuhs of BOl are compatible with a version of Keinigs' formula where a — r/tJ ■ B/B 

2 

is equal to —rjj ■ b/B ; see Eq. (48) of BOl. However, for the dynamo to be excited, the rjt term must always 
be slightly smaller than the a term. Therefore the sign implied for a is not affected by this generalization. 
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then a will be negative (Brandenburg 1998). The same result was also obtained by Riidiger 
& Pipin (2000) under the assumption that the turbulence is driven primarily by small scale 
magnetic fields. 

Next we compare quantitatively the results of the simulations with a simple mean field 
model. Here the negative sign of a follows directly from the fact that the dynamo wave is 
found to migrate away from the disc plane. 



3 Comparison with conventional mean field dynamos 

In mean-field dynamo theory one uses the averaged induction equation, 

dB 



dt 



V X {ux B + Su X 6B), (1) 



where overbars denote horizontal averages and the deltas denote fiuctuations. Our horizontal 
averages may be associated with ensemble averages in analytical studies. In astrophysical bodies 
with large scale shear, u = u{x), e.g. differential rotation, a mean component of the magnetic 
field is generated in the direction of the shear. However, this only works as long as there is 
a mechanism replenishing the magnetic field component in the cross stream direction. Such a 
field could arise from the correlation term 6u X 6B, the "turbulent" electromotive force, which 
we denote by £. Early work since the 1960s and 1970s has established the following form for 
£: _ _ 

Si = aijBj + r]ijkdBj/dxk (2) 

(Roberts & Soward 1975, Krause & Radler 1980, see also Steenbeck, Krause & Radler 1966 for 
an early reference, and Parker 1955 for the original formulation of the idea). In this representa- 
tion one assumes that the averaged magnetic field B is smooth so that higher derivatives in the 
Taylor expansion (^ do not enter. It is then assumed that the functional relationship between 
S and B is local. We return to the applicability of a local relationship in the next section. 
An apparently less general formulation of the diffusion term in equation (^ is 

Si = a.jBj - r]*pj, (3) 

where diffusion works only via the current density J = V x S (in units where the magnetic 
permeability is unity). However, the two cases are here equivalent, because in equation (^) 
only the index k = z gives a nonvanishing contribution (horizontal derivatives of horizontal 
averages vanish in our case!). Therefore the matrices rjijz and ri*j contain the same information. 
In particular. 



Vxx Vxy \ _ Vxyz 'V. 



Ixxz 



. ^ , , _ , (4) 

- Vyx Vyy ' \ Vyyz Vyxz / 

The advantage of formulation (|^) is that it is straightforward to ensure that turbulent diffusion 
does indeed lead to a decrease of magnetic energy. This is the case when the matrix r]*j is positive 
definite, i.e. rj*yJiJj > 0. A necessary condition for this is r^*^ > and rj*y > 0. Of particular 
interest is the case where 7]*^ is diagonal, i.e. 77*^^ and r]*^ are assumed to vanish. Below we shall 



7 



consider both representations, (H) and (|^) with rj*y = rj*^ = 0, and correspondingly we use 
different notation in the two and ri*j, respectively. We emphasize that equation (^) 

does not hold if the off-diagonal components of 7]*^ are put to zero, whilst rj^z is allowed to have 
all four components different from zero. 

In principle there could be an additional gradient term in the expression for S that would 
not affect the evolution of the field, but would modify the fitting procedure for the transport 
coefficient. (We are grateful to N. Kleeorin and I. Rogachevskii, who attracted our attention 
to this possibility.) In the present case where averages depend only on z, such gradient terms 
correspond simply to a constant vector, Eq] see, e.g., BD. However, in the present case the 
horizontal components of £ are antisymmetric about the midplane and have to vanish there. 
Since the terms on the right hand side of equations (|^) and are also zero on z = we must 
have Eq = 0. 

Now in order that the poloidal field can be replenished it is important to have a nonvanishing 
component of ayy. The magnitude of ayy has to be sufficiently large to overcome the effects of 
(turbulent) diffusion that arise from rjteijk, the isotropic part of rjijk- 

In BNST95, and later in Brandenburg & Donner 1997, the value of ayy was roughly esti- 
mated in the following way. It turned out that £y and By are correlated. The slope of the 
least-square fit gives a first estimate for ayy. It was found that ayy is negative in the upper disc 
plane and positive in the lower, and its magnitude was about 0.00ir2if, where Q = 27r/Trot is 
the local angular velocity and H the gaussian density scale height of the disc. In these papers 
no explicit estimate for rjt was given, but one would expect that rjt is comparable to the tur- 
bulent kinematic viscosity Ut, which was found to be approximately O.OOSfiif^. The negative 
sign of a in the upper disc plane was confirmed independently by Ziegler & Riidiger (2000), 
who used a different code. Their simulations also produced a large scale field. On the other 
hand, the simulations of Miller & Stone (2000) for a taller box did not show evidence for large 
scale dynamo action. It is possible that in their case the effective magnetic Reynolds number 
is larger and therefore the time scale after which a large scale field can be established would 
be longer. This effect would be even more pronounced if the numerical diffusion operator has 
hyper-diffusive properties (Brandenburg & Sarson 2002). 

In writing down the mean field equations for the horizontally averaged fields, we reiterate 
that because of periodic boundary conditions in the horizontal directions the averaged vertical 
field Bz, remains unchanged, and since Bz = initially, it remains so for all times. Thus, we 
have only two equations for the mean fields in the x and y directions, 

dB,__d£y 
dt ~ dz' 

dBy dS^ — 

_ = _ - ,SJB.. (6) 

The effect of molecular diffusion is here subsumed into the definition of £. On the boundaries 

we assume 

dB dB 

^ = ^ = on z = 0; B, = By = on z = Lz. (7) 
dz dz 
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This is equivalent to the conditions used in the three-dimensional simulations for the non- 
averaged magnetic field. The last term in equation arises from the local velocity shear, 
Uy{x) = —qflx. In the present case of keplerian rotation we have q = 3/2; for estimates of ayy 
for different values of q see Brandenburg & Donner (1997). The simplest parameterization that 
leads to dynamo action, balanced by diffusion, would be 

^. = r,.^, (8) 



£y = ayyBy-T]t^. (9) 



dz 

It is instructive to consider first free wave solutions of the form e^*"*"*^^, ignoring thus the 
boundaries. We now assume A; > 0. The dispersion relation is 

X = {l±i)\payyk\'/^-r]tk\ (10) 

where the two signs are respectively for positive and negative values of ayy. The marginally 
excited solution (ReA = 0) can be written as 



= Asmk{z - ct), By = V2A 



ayy 



sin[A;(z-ct) ± fTT], (11) 



where A is the amplitude (undetermined in linear theory), and c is the wave speed. 



_ ImA 

C = — = —ayy 



1/2 

= TVtk. (12) 



Note here that c > (as seen in the simulations, cf. Figure 0) requires ayy < 0. This is 
consistent with the sign of ayy obtained earlier by means of correlating Sy with By (BNST95, 
Brandenburg & Donner 1997). 

In Figure |]we present the result of a numerical integration of (|)-(||), where the boundaries 
are now taken into account. We chose 

ayy = -O.OOlOz, r/t = 0.005^i^^ (13) 

with Lz = 2H, and normalize to = = 1, so the rotational period is then Trot = 27r. 

The model reproduces roughly the average behavior of the magnetic field in the three- 
dimensional simulation. Of course, the simulation shows strong fiuctuations, which are not 
explicitly present in mean field theory. In the simulations fiuctuations are natural because 
the fiow is turbulent and horizontal averages are only an approximation to ensemble averages. 
Nevertheless the agreement between the mean-field model and the turbulence simulation is 
striking and even quantitative: the cycle period is 2 x 15 rotational periods, the migration 
speed is positive, although about three times too fast (c ~ +0.07QH), and the phase difference 
between B^ and By is about 37r/4 (Fig. |). Also the amplitude ratio between B^ and By is 
similar for the simulation and the model (about 1/30). Thus, the only major disagreement is 
in the migration speed, which is too fast in the model compared to the simulation. 

Our estimates of turbulent transport coefficients in equation (plSj) have been purely phe- 
nomenological. What is missing is an analysis of the validity of the expressions (|]) or (^). This 
will be attempted in the following section. 
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Figure 3: Horizontally averaged radial {B^) and toroidal (By) magnetic fields as a function 
of time (in rotational periods) and height, z, as obtained from the one-dimensional dynamo 
model. Dotted contours denote negative values. 

4 Fitting the functional form of 8 

In this section we attempt a direct determination of turbulent transport coefficients. We now 
allow alpha and turbulent diffusivity to be tensors whose components depend either on z (local 
approach) or on the vertical scale or wavenumbcr (nonlocal approach). In both approaches we 
also consider the case where 1]*^ is purely diagonal. 

Here and elsewhere primes on Bi and £i denote 2;-derivatives and dots denote time deriva- 
tives. We also drop the overbars on B^ and By, which denote now horizontally averaged fields. 
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Figure 4: Comparison of (dotted lines) and By (solid lines) as functions of time for z = 1 for 
the simulation and the model. Note the similarity in the phase shifts between and By in the 
two cases. In both panels B^ has been scaled up by a factor of 20. Note that both the model 
and the simulation have a similar amplitude ratio between the two fields. The horizontally 
averaged toroidal magnetic field in the simulation is given in units of the equipartition value, 
Beq = (47rpn^)^/^. The model on the other hand is linear and the amplitude is therefore 
undetermined. The field is growing in time, showing that the mean-field dynamo model is 
supercritical. 



4.1 A local formulation 

The underlying turbulence is anisotropic and therefore alpha-effect and turbulent diffusivity 
are really tensors. From @ we have 

£x = <^xxBx + axyBy + r]xxzB'^ + VxyzB'y, (14) 

~ ^yxBx + OLyyBy + TjyxzB^ + TjyyzBy. (15) 

Since we know S and B at different times we can determine the eight coefficients aij and rjij^ by 
forming moments with Bx-, By, B'^, and By. This gives 2x4 equations for the 8 unknowns, aij 
and rjij^. The reason why the third index on rjij^ is always z is because only the z-derivatives of 
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Figure 5: The symbols denote the resuhs for aij{z) as obtained by solving ([T6|) and the solid 
line represents a five point running mean. 



Bi are nonvanishing. This in turn is because the horizontally averaged fields are independent 
of X and y (see previous section). 

The full system of equations can be written in the form of two matrix equations, 



x,y, 



with the matrix 



M 



(BxBy) 

{ByBy) 

(B'yBy) 



{B.B',) 

(ByB'J 

{B'.B'J 
{B'yB'J 



{B.B'y) \ 

(ByB'y) 

(BLB'y) 
(B'yB'y)) 



(16) 



(17) 



/ {BxBx) 

(ByB^) 

{B'M 

\{B'yB,) 

which is the same for both equations {i = 1 and i = 2), and the vectors 



^(0 



The averages are taken over time. In Figures ^ and ^ we show the results respectively for the 
coefficients aij and rjij^ as functions of z. We also plot as a solid line a five point running mean 
of the data. 

Here and below our error bars are a measure of the stability of the fit when only part of the 
time series is used. In practice we divide the data set into two parts and calculate the transport 



/ {S^B,) ] 






i^iBy) 




(Xiy 


{£.B'J 




V {S^B'y) ) 




\ Viyz / 
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Figure 6: The symbols denote the results for rjij^^z) as obtained by solving (0) and the solid 
line represents a five point running mean. The two horizontal dash-dotted lines give the average 
values. 

coefficients separately for each subset. The error bars cover then the range of values obtained 
by using all or only part of the data. The error bars would probably shrink if we extended the 
data set, but some of the noise is due to fluctuations from cycle to cycle and therefore physical. 

Near the top boundary the data points deviate signiflcantly from the relatively smooth trend 
seen in the data away from the boundary. Therefore we consider the data near the boundary 
as uncertain and, when calculating averages and running means, we ignore all points within 
four mesh zones near the boundary. The vertically averaged values of the components of aij 
and Tjijz are 

_/ 0.010 -0.009 \ _/ -0.007 0.003 \ 2 /.q^ 

"'^■"U.ooi -o.ooj ' '''^■^ ~ I 0.000 o.oooy ^ ' 

Note that ayy = —O.OOlflH and rj^yz = O.OOSQH'^, in rough agreement with the results found 
earlier from other considerations. However, —rjyxz is at least ten times smaller and perhaps 
even negative. This coefficient is responsible for the diffusion of Bx-, i.e. B^ = ... — rjyxzB'J.. If 
—Vyxz < 0, this would indicate that some higher order terms (hopefully with the right sign!) 
would need to be restored for stabilization. Typically, an infinite series of further terms could 
then become important (Dittrich et al. 1984, Elperin et al. 2000), making this whole approach 
difficult to use in practice. We address this difficulty in §§ |4.3| and |4.4| by adopting a nonlocal 
approach. 
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Figure 7: The symbols denote the results for aij{z) as obtained by solving (|20|) and the solid 
line represents a five point running mean. Here and below the error bars indicate the stability 
of the fit when only a subset of the data is used. 



4.2 A local formulation using a diagonal diffusion tensor 

We now adopt the J-formulation for the diffusion term, i.e. equation (0), using a diagonal 
diffusion tensor, i.e. ri*y = t]*^ = 0. In that case we have to solve two systems of three (instead 
of four) equations. 



x,y, 



(20) 



where 

/ {BxBx) {BxBy) 

M(^) = {ByBx) {ByBy) 

\—{JiBx) —{JiBy) 

is now different for the two equations {i = 1 and i = 2). 
i is implied!) The vectors E^^^ are given by 



-{BxJ^)\ 

-{ByJi) (21) 
{JiJi) J 

(Here and below, no summation over 



(22) 



= {£,By) , c(^) = 

Again, the averages are taken over time. In Figures ^ and || we show the results respectively 
for the coefficients a^j and 7]*^ as functions of z. We also plot as a solid line a five point running 
mean of the data. 

Note that both rj*^ and 7]*^ are now positive on average (indicated by the dash-dotted line) . 
Nevertheless, ri*y, which is responsible for the diffusion of B^, shows still large fluctuations and 
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Figure 8: The symbols denote the results for r/j* (z) as obtained by solving ( pOD and the solid 
line represents a five point running mean. The average values are indicated by horizontal 
dash-dotted lines. 

can locally still be negative. Taken at face value, such excursions of r/*^ to negative values would 
cause a catastrophic growth of small-scale structures in a mean-field model. We note that the 
relation (^) between the different components of r/^*- and rjij^ still holds approximately, even 
though this relation is no longer strictly valid if, as in the present case, the fit to the diagonal 
components of 1]*^ is obtained under the restriction that the off-diagonal components vanish. 

4.3 A nonlocal formulation 

We recall that equation (H) corresponds really to a Taylor expansion of an underlying integral 
kernel in terms of derivatives of delta functions (see, e.g., Hasler, Kaisig & Riidiger 1995; see 
also Nicklaus & Stix 1988). Therefore we now write equations (0) and (|l^) as a convolution 
in the form 

*B^ + a^y *By + T],^,^^ * B'^ + r]^y^ * By, (23) 

£y = OiyX * B^+ ayy * By+ Tjy^^ * B'^+ Tjyy^ * By, (24) 

where aij{z, z' , t) and riijz{z, z', t) are now integral kernels and the asterisks refer to a convolu- 
tion, 

Oij * Bj = / aij{z, z')Bj{z' ,t)dz' , (25) 



T]ij, * = j^' rii,,{z, z')B'j{z', t)dz'. (26) 
We solve for the 8 kernels using Fourier transformation. We introduce 

e' = / g' ( z) cos kzdz, (27) 
Jo 

Bi = / Bi{z) cos kzdz, (28) 
Jo 
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where 



k = kn = {n + l/2)Ti/L„ n = 0,l, ... (29) 

Note that we use here the Fourier transform of the derivative of £, which is more easily obtained 
from the simulations. We assume that z') and fjijziz, z') have the form 



where ko = n/ {2Lz 



We also need 



sin koz {2/Lz) ^sin kz sinkz' aij{k), 

k 

{2/L;,) ^sin/czsin kz' fjijz{k), 

k 



Hk) 



Bi{z) sin kfjZ sin kz dz, 



and the corresponding inverse transform in terms of Bi{k), 

Bi{z) = {2/ L,){smkQz)-^Y.^mkzBi{k). 



(30) 
(31) 

(32) 
(33) 



The sin k^z term in (^) ensures that the a-effect is antisymmetric about the equatorial plane. 
For example, if C(ij{k) is independent of k, then the integral in equation (^) simply becomes 



a multiplication, so aij * Bj = sin fcoz«;>,(z), where alj' = aij{k) is a constant, independent 



(0) 



of k. We are now able to write (^) and ( p^ in the following form 

kOixxBx ~\~ k&xyBy k fjxxzBx k flxyzB 

Sy kCXyxBx ~t~ kCXyyBy /C ^jyXzBx k fjyyzBy 

Again, we can solve for the 8 coefficients dij and fjijz by computing moments (separately 
for each value of k). We obtain two matrix equations of the form 



(34) 
(35) 



E^'\k) = M{k)C^'\k), t = x,y 



with the matrix 



M 



( {BxBx) (BxBy) ^ 

{ByBx) (ByBy) {ByBx) {ByBy) 

(BxBx) {BxBy) (BxBx) {BxBy) 
V (ByBx) {ByBy) {ByBx) {ByBy) J 

which is the same for both equations {i = 1 and i = 2), and the vectors 



(36) 



(37) 



/ i^lBx) \ 

{£lBy) 

is'A) 

V {£lBy) ) 



/ kdix \ 
kd. 



k 'f]ix2 
V k'^fjiyz J 



(38) 
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Figure 9: The symbols denote the results for aij{k) as obtained from (|36| ) and the solid hne 
represents a five point running mean ignoring the point k = ko. 



The averages are taken over time. In Figures ^ and ^ we show the results respectively for the 
coefficients aij and fjijz as functions of k. We also plot as a solid line a five point running mean 
of the data. On the boundaries the lines go through the original data points, however. 

In the traditional afi-dynamo the {y, ?/)-component is the most important one. However, 
diyy{k) shows substantial variations in fc, indicating possibly large uncertainties, although cxyy 
is mostly negative and the value for ki is around —0.001, similar to the value we used in 
equation For the diffusion coefficient fj^yz the value for /c = fci is positive, i.e. stabilizing. 
Furthermore, the value for k = ki is much larger than for the next higher values of k. One 
possible suggestion would be to assume that r]ijz{k) was proportional to However, this 

would correspond to diffusion of the algebraic form Bi = ... — Bi/T^^\ i.e. without any differential 
(or integral) operator. Here the superscript (i) on r indicates that the decay times are in general 
different for the x and y directions. Let us recall, however, that because the eigenvalues of the 
diffusion tensor are not all positive our fitting is still not self-consistent. 

4.4 A nonlocal formulation using a diagonal diffusion tensor 

Like in § [4.2| we now adopt the formulation using a diagonal diffusion tensor. We write 

^'x = ka^xBx + ka^yBy - r)*^ J^, (39) 

S'y = k&y^B^ + k&yyBy " f^lj'y, (40) 
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Figure 10: The symbols denote the results for fjij^ik) as obtained from (|36|) and the solid 
line represents a five point running mean. Note that the data points for the smallest value of 
k, k = ko, deviate strongly from those for ki and lie outside the plot (57xxz(^o) = O.ISQH^, 
Vxyziko) = -0m07nH^, fiyUko) = OmnH^, and fj*(ko) = 0.002ni/2)_ 



where = +k'^By and J'y = —k'^B^. Instead of equation ( p6| ) we have 



E^'\k) = M%)&'\k), t = x,y 



(41) 



with the matrix 



(ByBx) 



(BxBy) 
{ByBy) 

-{J[By) 



-{BJi)\ 

-{ByJD 



(42) 



(Again, here and below no summation over i is implied!) The vectors are 



/ {^IB^) 

{£lBy) 



k&i:, 

k&iy 



(43) 



In Figures |ri| and we show the results respectively for the coefficients Uij and r)*^ as functions 
of k. We also plot as a solid line a five point running mean of the data. 



Now both fil^ik) and V^yik) are positive in the range 1 ^ ^ 15. Like in §^]^ the restriction 



to a purely diagonal magnetic diffusion tensor has helped to make turbulent diffusion positive 
definite. For larger values of k, corresponding to smaller scales, the data are probably no 
longer reliable. For the smallest value of k, k = ko, the diffusion coefficient of By is negative. 
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Figure 11: The symbols denote the results for aij{k) as obtained from (^) and the solid line 
represents a five point running mean ignoring the point k = k^. 
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Figure 12: The symbols denote the results for riij{k) as obtained from ( |^ ) and the solid line 
represents a five point running mean. Note that the data points for the smallest value of fc, 
k = ko, deviate strongly from those for ki and lie outside the plot {fjxxi^o) = —O.OOiQH^ and 



corresponding to a destabilization of the field on large scales. While this is of course what 
is seen in the simulation, it is unclear whether this is really due to a negative diffusion on 
large scales (like in the Kuramoto-Sivashinsky equation with negative diffusion and positive 
hyperdiffusion) , or just a natural manifestation of the a-effect or an inverse cascade. Let us 
note that the negative diffusion in the largest scale does not lead to any catastrophic behavior 
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of the solution. For recent work on the negative diffusion effect see Zhehgovsky et al. (2001). 



4.5 Comparative remarks 

In the previous subsections we have seen that a-effect and turbulent diffusion are not only always 
tensors, but their components show also systematic spatial variations and, more importantly, 
they are dominated by only the smallest wavenumbers. It does not seem feasible, however, to 
determine simultaneously spatial variations and wavenumber dependence. We emphasize that 
we do not advocate that the model of § |4.4| is better than that of § [4.2| or even for example. 



Instead, elements of both approaches should preferentially be taken into consideration. 

Generally, we do find, however, that the restriction to only diagonal components of the 
r}*j{z) or fj*-{k) tensors is to be preferred. This is because, on physical grounds, these diagonal 



components should be positive. If one does allow for off-diagonal components, as in 

and ^73| , one finds (as expected) somewhat different results {—rjyxz is different from ri*y, for 

example) and diffusion of Br^ field is then mostly negative. 

We wish to emphasize the importance of simultaneously determining the various transport 
coefficients, as we have attempted here. To our knowledge, the only other time this was done 
was in BOl where scalar a and r]t quenchings have been determined. In BOl, both direct 
determination and fitting to a model were used to obtain transport coefficients. Of course, 
fitting to a model is a much more stable and reliable procedure, but it is necessarily model 
dependent and one still has to assess whether or not the model is actually consistent with the 
data. 

An important complication arises from the fact that the actual electromotive force contains a 
strongly fluctuating component which is physical. This affects the results obtained from a direct 
determination of transport coefficients and, consequently, many of the detailed features cannot 
be physically meaningful. To our knowledge, dynamo models with nonlocal (wavenumber 
dependent) transport coefficients were never considered before. In order to have a preliminary 
assessment of such models, and to obtain insight as to which components are crucial, we now 
consider simple model calculations with wavenumber dependent transport coefficients. 



5 Model calculations 

The direct determination of nonlocal transport coefficients used in §^ suggests that nonlocal 
effects could be important. In fact, the magnitudes of dij and r}*^ tend to decrease with k. A 
A;~^ dependence of fj*^ and fi*y would cancel the diffusion operator and would yield a simple 
(and again local!) — l/r damping term. The problem with this type of damping is clear when 
looking at the dispersion relation (|TD|). For simplicity we assume a scalar magnetic diffusivity, 
Vlxi^) = Vyyi^) = Vti^)^ ^^"^ th^t Vti^) = (tP)'^ we have 

X{k) = {l±z)\lQa,y{k)k\y^-r-\ (44) 

so for &yy{k) = const, we would have A(A;) > for some k > kcru- In order to prevent such an 
'ultraviolet catastrophe' we have to require that ayy{k) —>■ for large values of k. On the other 
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Figure 13: Model with ayy{k) = -0.001(1 - k/k^^^), k < kyaa^ = 4, r^^ = 0.01. 

hand, the product ayy{k)k should reach its maximum for relatively small values of k, because 
the simulations (Fig. |l]) suggested no systematic high-/c dependence. 

Looking at Figure |TT|, Ciyy{k) fluctuates essentially about zero, although it does show a 
negative peak at = 37r/4 k, 2.4. Although \axx\ is much larger than it is relatively 

unimportant because it multiplies B^-, which is small, and a^x only contributes to generating 
By field, which is more effectively generated by the shear term. Given that it is only ayy{k) 
which contributes to aQ-tjpe dynamo action, we need to model this term somehow. We do 
this by assuming a linear dependence on k for A; < 4 of the form 

ayy{k) = -Omi{l-k/k^.,^) for A;</e^ax = 4. (45) 
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For ?7t(fc) we assume 

fi,{k) = (rk^y^ with = 0.01. (46) 

The resulting butterfly diagram is quite reasonable; see Fig. |r^. The main difference compared 
with the local model shown in Figure ^ is that the dynamo wave migrates now slower. The fast 
migration was the main shortcoming of the local model, where field reversals happened almost 
simultaneously both at the midplane and at the upper boundary. 

The assumption of a scalar fit{k) is made purely for simplicity. For comparison we have also 
considered a model with Vyy{k) = 0.1fi*^{k) (as suggested by Figure 0) and T)l.^{k) = {rk"^)"^ 
with = 0.01 (as before). In order that the model is still only marginally excited, we had to 
reduce &yy{k) by a uniform scaling factor of 5. Except for an increase of the cycle period by 
about a factor of 2, the resulting butterfly diagram was almost the same as in (Fig. [1^). 

Given that (i) this approximation reproduces phenomenology, (ii) is based on a fit to the 
actual electromotive force (see Figures |Tl| and [1^), and (iii) the dissipation is positive definite, 
we consider equations (^) and (|46|) as the most reasonable parameterization if one wanted to 
use mean-field theory. 

The results for many of the components of the a- and r^-tensors are rather noisy. However, 
the fact that the obtained values for Oyy, ayy{k), and rjt and (rA;^)~^, are all of the same sign 
and of comparable magnitude suggests this is a stable result of the analysis. This result is also 
in agreement with what is required for a mean-field model to work. Firstly, the sign of cxyy has 
to be negative to make the dynamo wave move with a positive migration speed. Secondly, the 
magnitude of ayy has to be around 10~^ so that the dynamo period is around 30 rotational 
periods. Finally, rji has to be around 10~^ so that the dynamo is just weakly supercritical. 



6 Noise in the transport coefficients 

The plots for both local and nonlocal approaches show strong fluctuations in the various trans- 
port coefficients. This is caused by strong noise, which is much more pronounced in the electro- 
motive force than in the resulting mean field. It would therefore be surprising if the a-effect was 
really responsible for generating the much less noisy large scale magnetic field in the simulation. 
However, it is important to note that even a completely noisy a-effect with no net component 
can amplify magnetic fields (Moffatt 1978). Furthermore, in the presence of large scale shear 
it is possible to generate also a large scale field (Vishniac & Brandenburg 1997). 

Given the high level of noise found in the transport coefficients it is useful to elaborate 
further on this point. In the following we present a model calculation where the effect of strong 
noise is included. We solve equations (|D-(|D and replace equation (|13D by 

ayy = -O.OOmz + anN, r]t = 0.0055fi//^ (47) 

where on is a coefficient and N{z, t) is gaussian noise in z and t, normalized to a root-mean- 
square value of unity. The other components of aij and rj*^ are ignored because we want to 
examine the effect of noise and take therefore the simplest possible model that was presented 



already in Figure Q. In Figure |14| we present the result for = 0.01, corresponding to a 
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Figure 14: Mean-field calculation with = 0.01, ao = —0.001, and rjo = 0.005. 



noise level that exceeds the mean alpha-effect by a factor ten or more, especially close to the 
midplane, where the coherent alpha-effect vanishes. 

The resulting field structure looks rather similar to the mean field obtained by averaging 
the data from the original three-dimensional simulations. The noise level exceeds the mean 
(coherent) alpha-effect by more than a factor of ten, which is comparable to or even larger than 
what has been suggested in connection with the solar dynamo (Choudhuri 1992, Moss et al. 
1992, Hoyng et al. 1994, Otmianowska-Mazur et al. 1997). The main point we want to make 
here is that a high noise level in the mean- field transport coefficients (§^) is quite natural and 
not simply an indication of inaccurate measurements. 
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7 Comments on nonlinear quenching 



The issue of magnetic quenching of turbulent transport coefficients has not been addressed in 
the present paper. One must expect that all components of the a and r] tensors depend on 
magnetic field strength and field structure. In principle our approach could be modified to allow 
for such nonlinearities. The usual assumption is that, if the turbulent transport coefficients are 
caused by turbulent convection, for example, a and rjt would decrease with increasing field 
strength. On the other hand, if the turbulent transport coefficients are caused by flow fields 
which themselves are magnetically driven, for example by magnetic (e.g., Balbus-Hawley and 
Parker) instabilities, then a and ?7t may actually increase with increasing field strength (see 
Brandenburg 1997 for a recent discussion on that). There is even some observational evidence 
for this somewhat unusual proposal (Brandenburg, Saar & Turpin 1998, Saar & Brandenburg 
1999). Thus, in addition to allowing a and r]t to be scale dependent, one should also allow 
for the possibility that a and r]t may actually increase with magnetic field strength. However, 
these type of behaviors must still be compatible with helicity conservation. 

In the case of a fully periodic domain the situation is now fairly well understood. For exam- 
ple, we know that the correct helicity limited growth of the large scale field can be reproduced 
if both a and rjt are equally strongly quenched by the mean magnetic field, and the quench- 
ing is of the form -\- asB /B^^), where B^q is the equipartition field strength, and the 
coefficient as is proportional to the magnetic Reynolds number. This type of 'catastrophic' 
quenching goes back to the early paper by Vainshtein & Cattaneo (1992), and was confirmed 
using simulations (Cattaneo & Hughes 1996), but it is now clear that this behavior is primarily a 
consequence of magnetic helicity conservation (Blackman & Field 2000, BOl), and independent 
of the type of nonhnearity (Brandenburg & Subramanian 2000). Evidence for 'catastrophic' rjt- 
quenching was first found by Cattaneo & Vainshtein (1991) using two-dimensional simulations. 
Three-dimensional simulations of isotropic helical turbulence in a box yielded a resistively slow 
saturation behavior of a which is well described by a fit of the form B ~ 1 — exp 
where rj is the microscopic magnetic diffusivity, k is the wavenumber of the large scale field, 
and At is the time after which the small scale field has saturated. The necessity to obey this 
fit is sufficiently stringent to exclude for example cubic quenching behavior (e.g., Moffatt 1972, 

2 

Riidiger 1974, Riidiger & Kitchatinov 1993), or quenching of the form 1 — obB / B^^. As was 
recently pointed out by Field & Blackman (2002), a type of quenching that satisfies the mag- 
netic helicity equation exactly and that is therefore also compatible with resistively limited 
growth is dynamical quenching (Kleeorin, Rogachevskii & Ruzmaikin 1995). According to this 
description, a and r]i are no longer catastrophically quenched. In other words, the earlier con- 
clusion of BOl is not compulsory because of a degeneracy of the dynamo equations in the fully 
helical case. 

In a nonperiodic domain the situation is quite different, because there is then the possibility 
of a magnetic helicity flux out of the domain. This magnetic helicity loss is associated with 
a loss of magnetic energy. In the model of BD much of the field that was lost through the 
boundaries was of large scale, which made the large scale dynamo process less efficient. This 
corresponds to a less drastic quenching, because is now only proportional to the square root 
of the magnetic Reynolds number (BD). 
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The presence of shear affects the saturation results in a different way. Field amplification 
by shear is quite independent of the magnetic helicity effect, but this concerns only the toroidal 
field amplification. For the dynamo to work one still needs to generate poloidal field in order 
to make use of the shear. Thus, saturation of the large scale field still occurs on a resistive 
time scale, but larger field amplitudes are now possible (Brandenburg, Bigazzi & Subramanian 
2001) Because of shear, the field evolves in an oscillatory fashion (as predicted by mean-field 
af2-dynamo theory; e.g. Moffatt 1978). It is at present not entirely clear, however, whether 
the time scale for this oscillation is resistive or dynamical. 



8 Conclusions 

In the present paper we have demonstrated that the results of the three-dimensional numerical 
nonlinear dynamo simulations of BNST95 can be fitted by a mean-field dynamo model with 
properly chosen parameterizations of transport coefficients. The parameterization found from 
such a fitting procedure is far from that of naive kinematic mean-field dynamo theory. In 
fact, our investigations point toward the possibility that turbulent transport coefficients might 
be wavenumber dependent. This possibility was first explored by Moffatt (1983), who used 
a renormalization group approach to obtain differential equations for a and r^t as functions 
of wavenumber. We find that the strongest contribution to a and r^t tends to come from the 
largest scales in the system. Unfortunately the scatter in the data is considerable and only 
data for the first few wavenumbers seem to be significant. However, simple representations of 
the form ( ^5] ) and ( |46| ) seem to reproduce the original data reasonably well (compare Fig. |13 



with Fig. |T]). These representations are consistent with fits to the electromotive force and, more 
importantly, the corresponding energy dissipation is now positive. 

Note that we have restricted the diffusion tensor that multiplies the current to be diagonal. 
This restriction means that the turbulent diffusion does not mix x and y components of the 
mean magnetic field. Such kind of 'cross-talk' is absent in dynamo models with short-correlated 
random velocity fields (Molchanov et al. 1985), although it may exist in other models (Radler's 
f2 X J-effect is an example; see Krause & Radler 1980). We note that our estimates for the 
a-tensor are not significantly modified when the diffusion tensor is restricted to be diagonal. 

Comparing the results of §§^ and 0!-|0, several similar aspects can be recovered in each 



case. Firstly, oiyy is negative in the upper hemisphere (or disc plane) and its magnitude, in 
natural units, is around 10"^ at z ^ H (see Figures ^ and |^ or at small wave numbers 
{k = ki), i.e. at large scales (see Figures ^ and 0). Secondly, the turbulent magnetic diffusivity 
in the streamwise direction (y) is such that the corresponding mean-field dynamo is marginally 
excited. The turbulent magnetic diffusivity in the cross-stream direction (x) is however much 
smaller than in the streamwise direction (cf. the two panels in Figure |^). This is again reflected 
in the nonlocal approach (cf. the two panels in Figure |12D, but here the diffusion coefficient of 
Bx is more noisy and, at the minimum wavenumber ko, the diffusion coefficient of becomes 
comparable to that of By, but it stays at least positive. We regard similarities in the results 
of the various approaches as an indication of their robustness (e.g., the values of ayy and rj*^). 
On the other hand, dissimilarities and noisy results, as well as large error bars, indicate that 
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results are not reliable (e.g. rjyy and rjyy)- Positivity of the diffusion coefficients is an additional 
criterion for reliability. 

We should stress that we expect realistic transport coefficients to show both height depen- 
dence as well as scale dependence. In addition, the a-effect is a dynamical one that satisfies an 
explicitly time-dependent equation (Rogachevskii & Kleeorin 2001). It is however rather diffi- 
cult to determine all these aspects simultaneously, which is why we have looked at each of them 
separately. Indeed, our goal is therefore not to distinguish between the various approaches, but 
rather to present a first assessment of the presence of each of these effects (height dependence, 
scale dependence) , all of which should be present simultaneously. 

We stress that our final model is based on very noisy fits for Uyy. We feel that it is im- 
portant to acknowledge the reality and physical significance of unsteady and random transport 
coefficients and to investigate models with stochastic alpha. Note also that there arc some 
models showing the possibility of the generation of large scale magnetic patterns in a flow with 
a random (incoherent) a-effect (Vishniac & Brandenburg 1997). However, in order to reproduce 
the butterfly (space-time) diagram obtained from the three-dimensional simulations a coherent 
(non- noisy) a-effect must still be present. The present results suggest that the simulations are 
reproduced with a coherent a-effect that is only about 10% of the incoherent a-effect in parts 
of the domain. 

The idea that only the smallest wavenumbers contribute to the transport coefficients may 
well prove to be a reasonable representation for astrophysical dynamos. Already now there 
is some evidence for this proposal: so far only models that include the lowest wavenumber 
are able to reproduce stellar cycle frequencies that decrease sufficiently rapidly with increasing 
rotational frequency (Brandenburg, Saar & Turpin 1998). By contrast, models with transport 
coefficients that are independent of wavenumber (Tobias 1998) show that the cycle frequency 
decreases too slowly with increasing rotational frequency. 

Although the present work was motivated by astrophysical and geophysical applications, 
the analysis presented here may just as well apply to homogeneous dynamos in the laboratory 
such as the Karlsruhe experiment (Stieglitz & Miiller 2001). The results from these experiments 
can now accurately be tested against simulations (Tilgner 2000) and mean-field theory (Radler 
et al. 1998). Due to the presence of boundaries, however, the a and r] tensors have to be 
replaced by integral kernels (Stcfani, Gcrbcth & Radler 2000). One may therefore expect that 
the approach explored in the present paper to determine integral kernels for the a and r] tensors 
will soon gain in importance. 
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